Dark Matter Annihilation Signatures from Electroweak Bremsstrahlung 
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We examine observational signatures of dark matter (DM) annihilation in the Milky Way aris- 
ing from electroweak bremsstrahlung contributions to the annihilation cross section. It has been 
known for some time that photon bremsstrahlung may significantly boost DM annihilation yields. 
Recently, we have shown that electroweak bremsstrahlung of W and Z gauge bosons can be the 
dominant annihilation channel in some popular models with helicity-suppressed 2^2 annihilation. 
IV/Z-bremsstrahlung is particularly interesting because the gauge bosons produced via annihilation 
subsequently decay to produce large correlated fluxes of electrons, positrons, neutrinos, hadrons (in- 
cluding antiprotons) and gamma rays, which are all of importance in indirect DM searches. Here we 
calculate the spectra of stable annihilation products produced via 7/ bremsstrahlung. After 
modifying the fluxes to account for the propagation through the Galaxy, we set upper bounds on 
the annihilation cross section via a comparison with observational data. We show that stringent 
cosmic ray antiproton limits preclude a sizable DM contribution to observed cosmic ray positron 
fluxes in the class of models for which the bremsstrahlung processes dominate. 

PACS numbers: 95.35.-|-d, 12.15.Lk, 95.85.Ry 



I. INTRODUCTION 

It has been firmly established that a significant frac- 
tion of the energy density of the Universe resides in the 
form of dark matter (DM), whose total abundance and 
approximate distribution has been inferred via its grav- 
itational infiuence [iHl]- However, we have as yet no 
direct detection, and the fundamental particle properties 
of DM remain unknown. A promising technique to ex- 
plore the particle nature of DM is to search for signatures 
of DM annihilation or decay. This is achieved by exam- 
ining fluxes of particles emanating from regions of DM 
concentration either in our own Galaxy or throughout 
the cosmos. 

Indirect DM detection has been the subject of much 
recent attention, due to measured cosmic ray excesses 
of positrons and electrons above those expected from 
conventional astrophysical processes. PAMELA has ob- 
served a sharp excess in the e+/(e^ + e^) fraction at en- 
ergies beyond approximately 10 GeV 4], without a cor- 
responding excess in the antiproton/proton data 0, Q, 
while Fermi and HESS have reported more modest ex- 
cesses in the (e~ -I- e"*") flux at energies of order 1 
TeV 0, Hi- These signals have led to a reexamination 
of positron production in nearby pulsars [To| . emis- 
sion from supernova remnants [11|, acceleration of e~^e~ 
in cosmic ray sources |12|^and propagation in conven- 
tional cosmic ray models [lj|. Alternatively, it has been 
proposed the excess electrons and positrons are not due 
to conventional astrophysics process, but arise instead 
from DM annihilation or decay in the Galactic halo. A 
plethora of DM models have been designed with this goal 
in mind (for a review see for example [ij] , and references 
therein). 



A viable resolution of the cosmic ray data by means 
of DM annihilation requires a large branching ratio to 
leptons. A large branching ratio to hadrons would make 
a contribution to cosmic ray antiproton fluxes, for which 
stringent observational bounds exist. Therefore, so called 
leptophilic models are preferred, in which DM couples (at 
tree level) only to leptons. However, for many scenarios 
in which the DM particle is a Majorana fermion, annihi- 
lation to light fermions is helicity suppressed (oc m^/s) 
in the s-wave contribution, and of course velocity sup- 
pressed (oc vj^j^.j) in the p-wave contribution. This is 
the case for popular DM candidates such as the neu- 
tralino of supersymmetric models, if ^/s ^ 2Mdm is be- 
low the W^W^ threshold; or Bino models with highly 
suppressed annihilation to W~^W~ and ZZ final states. 
Large boost factors would be required for such a scenario 
to explain any observed positron excess. 

It has long been known that pho ton bremsstrahlung 
can lift helicity suppressions |15l - [l7| . We have recently 
shown that helicity suppressions are also lifted by the 
bremsstrahlung of a or Z gauge boson [ll,[i^ (see also 
[20}). (Electroweak radiative corrections to DM annihila- 
tion, but without the observation of helicity unsuppres- 
sion, have been discussed in Ref. f2ll - [28| .) Where a helic- 
ity suppression is lifted, the cross sections for xx ~^ ^^1: 
XX ~^ ^J^d XX ~^ can all greatly exceed that for 
the lowest order process xX ~^ The bremsstrahlung 
processes thus may allow for the indirect detection of 
many DM models which would otherwise be helicity sup- 
pressed. 

Importantly, the decay of the W and Z gauge bosons 
inevitably leads to the production of secondary annihila- 
tion products, including gamma rays, hadrons, charged 
leptons and neutrinos, allowing multimessenger searches. 
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Note particularly, that even for DM models designed to 
be leptophilic, production of hadrons is unavoidable. (In 
fact, even for models in which on-shell production of W or 
Z gauge bosons is kinematically forbidden, some minimal 
hadron production is inescapable, due to loop processes, 
or the exchange of off-shell W or Z bosons.) 

In this paper we examine a simple example of a model 
which has a helicity-suppressed 2 — >■ 2 cross section, pro- 
vided by Ref. [l^ls^l- We shall show that in this model, 
the cross section required to produce positrons in suf- 
ficient quantity to account for the observed excess will 
lead to overproduction of antiprotons and gamma rays, 
and as such is ruled out as an explanation of the observed 
positron anomalies. Though our calculation is performed 
for the reference scenario of Ref. [1^ [sO] , we expect this 
conclusion to hold for all scenarios in which the domi- 
nant positron contributions arise from the 3-body W / Z- 
bremsstrahlung final states. 

We calculate the spectra of both primary and 
secondary particles from unsuppressed electroweak- 
bremsstrahlung annihilation processes, and calculate the 
expected spectra and fluxes at Earth for a given annihi- 
lation cross section. We compare the Earthly fluxes with 
observational data to determine an upper limit on the 
annihilation cross section. While our analysis techniques 
are conservative, there are large astrophysical uncertain- 
ties in the propagation of charged particles through galac- 
tic magnetic fields, and in the DM density profile which 
probably contains substructure. A rigorous treatment 
of these effects is beyond the scope of this work. Conse- 
quently, our constraints are illustrative of the upper limit 
on the cross section, but not robust. 



II. MODEL 

The example model we investigate is the Majorana DM 
version of the leptophilic model proposed in [2^. Here 
the DM consists of a gauge-singlet Majorana fermion x 
which annihilates to leptons via the interaction term 



X + h.c. = f{iyLv" - ^L^^)X + h.- 



(1) 

where / is a coupling constant, e is the S'?7(2)-invariant 
antisymmetric matrix, and (77"*", rf) form the new SU{2) 
doublet scalar which mediates the annihilation. For sim- 
plicity, we consider a coupling to the first generation of 
leptons only, and set / = for coupling to the (vg fJ.^)L 
and {vtT )l doublets. As described in [IMISIj the p- 
wave contribution to the lowest order annihilation pro- 
cess XX e+e^ is suppressed by ^ 10^®, while the 
s-wave contribution is proportional to [mi/M-^Y . This 
cross section is given by 
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where m; ~ and M^± = M^o have been assumed, and 
fi = M^/M^. The helicity-suppressed s-wave term is 
absent in the m; = limit, leaving only the u^-suppressed 
p-wave term. 

While it is well known that photon bremsstrahlung 
XX ~^ e"'"e~7 can lift this suppression (TB - flTl . 
Refs. [H, [1^ have shown that this is also the case 
for the electroweak bremsstrahlung channels XX ~^ 
e+e-Z, v^v^Z, e+i^eW', e-VeW+. For both W/Z and 
7 bremsstrahlung, the effect is most significant where 
the DM mass is nearly degenerate with the mass of the 
boson which mediates the annihilation process. In the 
high energy limit where the W/Z masses are negligi- 
ble, the cross sections for W/Z bremsstrahlung reduce 
to that for 7 bremsstrahlung, modulo different coupling 
constants. However, the respective sizes of the electro- 
magnetic and electroweak coupling constants imply that 
the ly/Z-strahlung cross section is a factor of several 
larger than that for 7-strahlung, 
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and thus 



V (Tw/Z-strahXung = 6.16w(Te 



(3) 
(4) 

(5) 
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At lower energies, the phase space for W/Z 
bremsstrahlung is somewhat reduced due to the 
effects of the finite W/Z masses. 

The bremsstrahlung cross section dominates over that 
for the lowest order 2 — )■ 2 process, provided that 
Mn does not greatly exceed M^. The ratio Rw = 
v<Je+v^w- /{"^^e+e-) was plotted in Ref. [l^, and is 
largest for ^ = [M^/M^Y = 1. However, the W- 
strahlung process dominates over the 2 — >■ 2 annihilation 
even if a mild hierarchy between M^^ and Af^ is assumed, 
with R\Y > 1 for fj, < 10. Therefore, when /i < 10, pro- 
duction of monoenergetic leptons will be subdominant to 
particles produced in the 3-body processes (both primary 
and through gauge boson decay), a feature which must 
be accounted for when analysing astrophysical signatures 
of these models. 

We use the cross sections calculated in [ 16l -[l8| in com- 
bination with the PYTHIA code [Sll, |32|| to determine 
the spectra of gamma rays, electrons, protons, and their 
antiparticles, per annihilation to the five 3-body final 
states listed above (one electromagnetic bremsstrahlung 
and four electroweak bremsstrahlung processes). After 
accounting for propagation effects for the charged parti- 
cles, we constrain these cross sections by comparing the 
observed flux with the calculated annihilation signal. 

Although the spectra of annihilation products which 
we show are unique to the particular model we have cho- 
sen, we expect the results to apply qualitatively to any 
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model where VF/Z-bremsstrahlung is the dominant an- 
nihilation mode (i.e. where helicity suppression of the 
2 — )■ 2 s-wave is lifted by electroweak bremsstrahlung) . 



III. CROSS SECTION CHANNEL 

As mentioned, the model we examine is lep- 
tophilic in the 2 — > 2 process, and therefore 
has five 3-body bremsstrahlung channels, XX ~^ 
e'^e~ Z, VeVeZ, e'^ VeW~ , e~ i>eW~^ , e"'"e~7, which simulta- 
neously contribute to the DM annihilation fluxes. The 
cross sections and spectra for the electroweak channels 
are specified in Ref. [l^, while those for the electro- 
magnetic channel are given in Refs. The total 
bremsstrahlung cross section is given by the sum 

+ va^+v^w- + WCre-p^iy+ + WCTg+e-^. (7) 

The total bremsstrahlung cross section is a factor of ^ 7.2 
larger than that for photon bremsstrahlung alone, due to 
the four W I Z channels, which are governed by somewhat 
larger coupling constants. 

We shall consider parameters for which the 
bremsstrahlung channels dominate the total cross 
section, so that I'O'Brcm — vatotai- We define 
the "branching ratio" for an individual channel 
i e {e^e~ Z,i>eVc.Z, e'^VeW~, e~i>eW~^, e~^e~j} as 



BRBrcni{i) 



VaBi: 



(8) 



The spectrum per annihilation, for any given annihilation 
product, k G {7, e~, e"*", i/, i?, p, p}, is then given by 



dNk 



dEk 



= ^ Si?Bicin(i) 



dNk 



dEk 



(9) 



Here, 



dNk 
dE,. 



is the spectrum per annihilation for 

per xX^i 

a given channel. The spectra for 7, e^, v and v include 
primary annihilation products and secondary annihila- 
tion products produced by gauge boson fragmentation. 
The spectra for p and p arise exclusively from gauge 
boson fragmentation. The branching ratios and spec- 
tra depend on the parameter ^ — {M^/M^Y . However, 
as long as the 3-body final states remain the dominant 
channel, the spectra (and thus the final results) have lit- 
tle dependence on this parameter. We show results for 
ji = {Mr^/My^i^ = 1.2, but results remain qualitatively 
unchanged when bremsstrahlung channels dominate over 
2 — T' 2 processes. 

Finally, the flux of a given annihilation product is 
schematically 



d4>k dNk 
dEk " 



(10) 



Brcm 



The detailed evaluation of these annihilation spectra is 
given below. 



IV. ANNIHILATION SPECTRA 

In order to place constraints on the cross section, 
we need the spectrum of stable particles {v,e~,p and 
their antiparticles, plus 7) produced per DM annihila- 
tion. As an example, we describe how we determine 
the spectrum of antiprotons per xX ~^ i'eS-'^W~ event; 
the technique is very similar for other secondary parti- 
cles, and other electroweak-bremsstrahlung annihilation 
channels. These partial spectra are then summed to form 
dNk 



dEk 

We use the PYTHIA code (3l|, |32| to find the spec- 
trum of antiprotons per W~ decay, dNp/ dE\^^^^^^. We 
produce a W~ boson in its rest frame by colliding an 
antimuon with a muon neutrino, with center of mass 
energy Mw, and turning off all processes other than 
^'^i'^{Vfj) -> W^. (Similarly, to produce the Z boson, 
we collide a e+e^ pair at CoM energy M^, leaving Z pro- 
duction as the only active process.) Unstable W~ decay 
products (mainly pions) themselves decay, finally leaving 
only neutrinos, electrons, protons and their antiparticles, 
plus gamma rays in the final state. These stable particles 
are placed in 2000 logarithmically-spaced energy bins. 
The final spectrum is found by averaging the PYTHIA 
spectra over 10,000 such events. 

We use this energy spectrum, in combination with 
the (or Z) energy distribution per annihilation, 

dNw/dj — {l/va){dva/dj), where 7 = Ew/Mw, to 
find the antiproton energy-spectrum per annihilation in 
the lab frame (see Appendix of [l^ for a derivation) 



dNp{E') _ 1 
dE' ~ 2 



d7 dN^ 



w 



d7 



^-dEdN, 
E P dE' ^''> 



with p = ^E^ -mj, (3-f ^ ^^7^ - 1, and E± = -fE' ± 
l3jp'. Or equivalcntly, 

dNp{E') _l r°° dE dNp r+ d-y dNw 



dE' 



p dE ^72 _ 1 d'y 



(12) 



with 7± = {EE' ±pp')/mp and p' = \J E'"^ — rUp. In the 

case of the gamma-ray, electron, positron, and neutrino 
spectra, we add the spectrum of primary annihilation 
products to the spectrum of secondaries from W decay, 
calculated as above. 

Figure [T] shows the total gamma-ray spectrum, as well 
as the relative contributions from primary and secondary 
annihilation products, clearly showing that the secondary 
gamma rays are subdominant, except at low energy. 
Figure [5] shows the same information for the positron 
spectra, including the relative contributions to the pri- 
mary positron spectrum from photon bremsstrahlung 
and electroweak bremsstrahlung, as well as the sec- 
ondary positrons produced via gauge boson decay. In 
Fig. 131 we show the total spectra per annihilation for 
electrons, neutrinos, protons, and gamma rays. Note 
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FIG. 1. Contributions to the gamma-ray spectrum per an- 
, from primary production in photon 



nmnation, 

dE 



bremsstrahlung (dotted, orange), and W/Z decay products 
(dashed, green), for = 300 GeV and M^/M^ = 1.2. The 
total gamma-ray spectrum is shown as a sohd curve (red). 
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FIG. 2. Contributions to the positron spectrum per annihi- 



lation, 



dE 



from primary production in electroweak 



bremsstrahlung channels (dotted, orange), primary produc- 
tion in the photon bremsstrahlung channel (dashed, green) 
and W/Z decay products (dot-dashed, blue), for — 300 
GeV and (Mn/M-^)'^ = 1.2. The total positron spectrum is 
shown as a solid curve (red). Note that the positron spectra 
from electroweak and photon bremsstrahlung have differing 
kinematic cutoffs due to the masses of the and Z bosons. 




0.01 



0.001 

0.4 0.6 

x=E/mx 

FIG. 3. Spectrum per annihilation of photons (solid, red), 
protons (dotted, orange), electrons (dashed, green) and neu- 
trinos (dot-dashed, blue), for = 300 GeV and M^/M^ = 
1.2. For protons, E is the kinetic energy. By CP-invariance, 
the particle and antiparticle spectra are the same, and an- 
tiparticles are not included in this figure. Note that the neu- 
trino spectrum includes primary electron neutrinos, and all 
flavors of secondary neutrinos. 



annihilation to a pair of gauge bosons. The photon spec- 
trum for the WIMP annihilation channel XX ~^ W^W~ , 
shown for example in Fig. 1 of Cembranos et al. [33| . 
has a somewhat softer gamma spectrum, while we have 
a somewhat harder spectrum with more higher-energy 
photons. This is to be expected, as gamma rays in 
Ref. [3^ arise only from decay of monoenergetic W 
bosons {Ew = M^), while the photon-bremsstrahlung 
process contributes a harder primary gamma-ray spec- 
trum. The reverse holds true for the proton spectra, 
as our 3-body M^/Z-bremsstrahlung process results in a 
broad distribution of W energies (the spectrum is shown 
in Fig. 3 of Ref. [3). In addition, the electron and 
neutrino spectra resulting from the xX ~^ W~^W~ pro- 
cess would be quite different to those for electroweak 
bremsstrahlung, given that for the latter it is the primary 
leptons (not the secondaries from W/Z fragmentation) 
that make the dominant contribution. Of course, the 
spectra of charged particles observed at Earth will differ 
from that at production, due to the effect of energy loss 
processes during propagation. We address these effects 
in the following section. 



that the electron/positron spectra from W/Z and pho- 
ton bremsstrahlung have differing kinematic cutoffs due 
to the masses of the W^ and Z bosons, leading to a kink 
near the endpoint in the electron/positron spectra. This 
feature is absent from the neutrino spectrum, as there 
is no contribution from photon bremsstrahlung. The 
neutrino spectrum includes contributions from primary 
electron neutrinos, and all flavors of secondary neutrinos 
from W^ /Z decay. (The flavor ratios of primary neutrino 
production are mo del- dependent, given by : : f^. 
We have assumed f^ = f^ = 0.) 

It is illuminating to compare our spectra to those for 



V. CONSTRAINTS 

In this section, we place conservative upper limits 
on the thermally-averaged self-annihilation cross section, 
(fcr) Brcm — (vct) total, for the Icptophilic model described 
above. We do this by following the same technique as 
in, for example, [13, [13, [s^- We compare the various 
predicted fluxes for a particular DM annihilation chan- 
nel with the relevant observational flux measurements. 
We make the conservative assumption that the entire ob- 
served flux comes from DM annihilation; in reality, as- 
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trophysical backgrounds are likely to contribute a large 
fraction of the observed fluxes. The upper limit on the 
cross section is then determined such that the DM anni- 
hilation does not exceed any of the observed fluxes. 

In calculating the constraints on (wcr)Brcmj we uti- 
lize the isotropic extragalactic gamma-ray flux measured 
by the Fermi collaboration |36| , the very-high-energy 
gamma-ray flux from H.E.S.S. |37|. the positron fraction 
from the PAMELA coUaboration the Fermi e+ -I- e~ 
flux , as well as the antiproton flux and antiproton-to- 
proton ratios updated by the PAMELA collaboration in 
Ref. . Wherever uncertainties in the flux are presented, 
we use the 1-a upper limit. Throughout, we use the 
commonly-adopted Navarro, Frenk and White (NEW) 
DM density profile [1^1, with local DM density given by 
Pq = 0.39 GeV /cm~^ [S^Q We properly account for 
proton and electron diffusion and energy loss as detailed 
below. 



A. Gamma Rays 

The isotropic diffuse gamma-ray flux will have contri- 
butions from both galactic and extragalactic DM annihi- 
lation. Although the galactic signal is expected to have 
a large directional dependence, there will be an under- 
lying isotropic component [43|. We include both contri- 
butions when computing constraints, though the galactic 
flux dominates over the extragalactic for the parameters 
of interest. In order to calculate the cosmic annihilation 
signal, we follow the technique set out in Refs. (il. |45|. 
while adopting the notation of {ssl, ^43!] . 

The isotropic gamma-ray flux from DM annihilations 
throughout the Universe is given by 

d^^ (va) c pI,, 



dE 



^ttHo Ml 
'^'^ f{z){l + zf dN^{E') 



h{z) 



dE' 



e-"("^^'dz, (13) 



where Hq = 70 km s ^ Mpc ^ is the Hubble parameter, 
Pav is the average DM density in the Universe, and h{z) = 
[{l + zf^DM+^hY''^- We assume Udm = 0.3, VIa = 0.7. 
The energy at production is E' , and the redshifted en- 
ergy IS E = E' l{\ -f z). We use the optical depth r(z, E) 
from , which accounts for attenuation of gamma rays 
as they propagate through the Universe. The factor /(z) 
accounts for the clustering of DM, which gives an en- 
hancement of the annihilation signal relative to a uni- 
verse in which matter was distributed homogeneously. 
Following psj . we parametrize the redshift dependence 
as logio(/(z)//o) = 0.9 [exp(-0.9z) - 1] - 0.16z. There 
is some debate as to the overall normalization of /(z). We 



P0 = 0.389±0.025 (l-cr) but could be a factor of ~ 2 higher 
or lower |40.--42] . 



adopt the NEW profile, which gives /o ~ 5 x 10^. This is 
a conservative choice for /(z), and inclusion of enhance- 
ments due to subhalos would only strengthen our results. 
See Ref. '46^ for a recent discussion of the clustering fac- 
tor. Our choice for /(z) closely corresponds to one of the 
smallest examples given there. 

For Galactic annihilation, we again follow the tech- 
nique of e.g. ^35i . i43i] . The flux of gamma rays per stera- 
dian from Galactic DM annihilation, in a direction at an 
angle i\} from the Galactic Center, is given by 



'Je 



1 (w) Ji'4j)dN^ 



2 47rM2 Jo 



dE 



(14) 



where 



= Jo 



P 



Rl - 2£i?,c cos iP + e]d£ (15) 



is the integral along the line of sight of the DM den- 
sity squared. Jq = l/[8.5kpc x (0.3 GeV cm~^)^] is an 
arbitrary normalization constant used to make J^iip) di- 
mensionless; it cancels from our final expression for the 
gamma-ray flux. Since we are calculating the isotropic 
signal, we use the minimum value for J^{ip), once again 
adopting the NEW profile. 

The Extragalactic Gamma-Ray Background (EGB) re- 
ported by Fermi in [36| is the isotropic component of 
the diffuse gamma-ray flux, with a number of potential 
contributing sources. It is obtained by subtracting the 
components of the gamma-ray flux with known origin 
from the total flux, observed away from the Galactic disk 
(Galactic latitude |&| > 10°). Hence, it is a flux likely 
to contain a contribution from either Galactic or extra- 
galactic DM annihilation. We compare our calculated 
isotropic signal, from both cosmic and Galactic annihi- 
lation, to this isotropic flux. We do this for each data 
energy bin, integrating the signal over the width of each 
bin in turn. 

We also compare the DM annihilation signal with 
H.E.S.S. observations of the very-high-energy gamma- 
ray flux from an angular region of 1° around the Galactic 
Center, excluding the Galactic Plane [13]. These observa- 
tions have both advantages and disadvantages when com- 
pared with the Fermi EGB observations. First, the two 
data sets cover different energy regimes. The H.E.S.S. 
data begin at around 300 GeV, and so are unable to 
constrain the cross section at lower DM masses. Con- 
versely, Fermi EGB data end at 100 GeV, such that con- 
straints on the cross section for larger DM masses are 
based on comparisons of the observed flux with the tail- 
end of the annihilation spectrum, leading to relatively 
weak constraints. In addition, the two data sets are 
from drastically different observation regions. The ex- 
pected isotropic DM signal is relatively small, leading to 
weaker constraints from the Fermi data; however, this is 
balanced by the small uncertainty between the various 
DM density profiles for this region. The opposite holds 
true for the H.E.S.S. observation region, where there is 
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a large expected DM signal, large background flux, and 
a very large uncertainty between competing density pro- 
files, since the observation region is small and close to 
the Galactic Center. While we report results using the 
NFW profile, these would be weakened by well over an 
order of magnitude if an isothermal profile were adopted. 

For these reasons it makes sense to constrain the an- 
nihilation rate using both data sets, emphasizing that 
while we find that the H.E.S.S. constraints are signifi- 
cantly stronger, they are subject to larger uncertainties. 
We follow the same technique as earlier, treating the flux 
from the Source region in Ref. [37| as an upper limit on 
the annihilation flux calculated using Eq. [T3] for each re- 
ported energy bin, allowing us to place an upper limit 
on {va) as a function of the DM mass. We again present 
results using the NFW profile, utilizing the average of 
J{tp) over the observation region as reported in Ref. [37| . 
Jav = 1604. 

Our resulting upper limits on {va) using both Fermi 
EGB and H.E.S.S. data are reported in Fig.H] 



B. Electrons and Positrons 

The flux of positrons (or electrons) at Earth from 
DM annihilation depends both on the propagation of the 
positrons through the turbulent galactic magnetic fields, 
and energy losses of the particles. One can solve the ap- 
plicable diffusion-energy loss equation to find a semiana- 
lytic form for the positron flux at Earth. We adopt the 
same notation and set of assumptions as in e.g. Ref [47| . 
which gives 



dE ~ SttMME) Je 



'PqV 



,dN, 



dE'^I{\D {E,E' j) 



(16) 

where pq is the local DM density, b{E) and / (Xd (E, E')) 
are the energy loss and "halo function" parameters re- 
spectively, and Xd{E,E') is the diffusion length be- 
tween the two energies E and E'. The energy loss 
of the positrons as they propagate through the Galac- 
tic medium is mainly due to synchrotron radiation and 
inverse Compton scattering, and is characterized by 
biE) « 10-i'^(i;/GeV)2 GeV s-\ The "halo function" 
/ (A^i {E, E')) is an astrophysical parametrization which 
encodes the dependence of the flux on the DM density 
profile, and on the model of positron diffusion due to 
Galactic magnetic fields. There is a degree of uncer- 
tainty in this function, and so we choose the 'medium' 
diffusion parameter set (while also showing results in the 
'min' and 'max' scenarios), and as usual, the NFW DM 
density profile. We use the numerical fit to / (Ad [E, E')) 
from [47[. 

Our signal is then compared with the total e^-f flux 
reported by the Fermi collaboration Q to find an upper 
limit on {va), by demanding that the signal integrated 
over the width of an energy bin be less than the total 



observed flux in that bin. 

We can also combine our positron flux with the Fermi 
data to find the positron fraction from DM annihila- 
tion. We compare this with the PAMELA data 4] for 
the positron fraction (/e+) to find an alternative upper 
limit on {va) by demanding 



e++e- 



</e 



(18) 



We compare the DM-related positron fraction with the 
observed PAMELA fraction in each of the four energy 
bins where the Fermi energy range overlaps the PAMELA 
energy range, integrating the DM-signal and observed 
Fermi fluxes over the width of the PAMELA energy bins. 
For this, we use the simple power-law fit to the Fermi 
data, valid between around 20 GeV and 1 TeV [l0| . 



dE 



= (175.40 ±6.09) X 10-"* (GeVcm^ ssr)-^ 
X (^/GeV)(-3-045±o.oo8) _ (19) 

Results are reported in Fig. |4l while results in the 'min- 
imum' and 'maximum' diffusion scenarios are shown in 
Fig. El 

C. Protons and Antiprotons 

The antiproton (or proton) fiux at Earth has a sim- 
ilar functional form to that for positrons, except that 
the energy losses for the antiprotons as they propagate 
to Earth are negligible, since rup ^ m^. (Note that al- 
though energy loss is negligible, diffusion is not.) Be- 
cause the energy for the antiprotons is the same as the 
injection energy, there is no need for an integration over 
energies E' at production. (Energy loss due to scatter- 
ing interactions or solar modulation are relevant only at 
low energies.) We again use the semianalytic function 
from [i^l to calculate the proton and antiproton signals 
at Earth from DM annihilation. 



d^p{K) _ {va)p%v 



dK 



R{K) 



(20) 



where K is the kinetic energy of the (anti)-proton, and 
R{K) is an astrophysics parametrization playing a sim- 
ilar role to I{\d {E,E')) from Sec. |VB1 Reference [13] 
provides a numerical fit to R{K) for several sets of prop- 
agation parameters, and we again use the 'medium' pa- 
rameter set, with the NFW DM density profile. 

We compare our antiproton flux with the total antipro- 
ton flux reported by the PAMELA collaboration [5|, en- 
ergy bin by energy bin. We can also constrain the cross 
section by demanding the ratio p/p due to antiprotons 
from DM annihilation not exceed the PAMELA p/p ra- 
tio 0] . This comparison requires the observed proton 
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flux. Following [13], we use the nucleon flux from (isj , 

^ll^ w 0.79 X 1.8(i;/GeV)-2 '^ (cm^sGeV)-!, where 
0.79 is the proton fraction of the total nucleon flux. We 
then simply demand ^p/^°^^ < fp/p, where /p/n is the 
PAMELA antiproton/proton flux ratio given in Q. En- 
ergy bins have been handled in the same way as the 
positron case, giving us the upper limit on (va) shown 
in Fig. m and results using the 'minimum' and 'maxi- 
mum' diffusion parameter sets are shown for comparison 
in Fig. [5l 



D. Neutrinos 

Following Ref . p3l | , we examine the neutrino flux aver- 
aged over the entire sky. The dominant contribution to 
the DM-signal will be from Galactic annihilations, and 
we ignore the subdominant cosmic annihilation signal. 
We make the approximation that 1/3 of all produced 
neutrinos will be observed as muon neutrinos. 



dE 3 dE ' ^ ' 

Note that the neutrino spectrum in Fig. [3] includes pri- 
mary electron neutrinos, and all flavors of secondary neu- 
trinos. We then compare the -I- signal from DM 
annihilation with the atmospheric v^, -\- flux, using the 
same technique as Ref. [43] . The calculation of the Galac- 
tic neutrino signal proceeds the same as for gamma rays, 
using Eq. [14] with the neutrino-annihilation spectrum in 
place of the gamma-ray spectrum. 

We average J{ip) over the entire sky, and find Jav — 5 
for the NFW profile. The signal is then compared with 
the background flux, integrating each over an energy bin 
width of AlogiQ E = 0.3. As expected, the resulting up- 
per limit on the cross section is significantly weaker than 
those calculated using the other annihilation products 
considered, and is only visible in the high-mass region of 
Fig. 21 Accordingly, the assumptions made in this sub- 
section concerning neutrino flavors are moot. 

As an alternative technique, one could calculate the 
flux of upward-going muons through the Earth induced 
by neutrinos from DM annihilation, and compare this 
with limits from the Super-Kamiokande (SuperK) exper- 
iment Reference f50] converts the SuperK limits 
into an upper limit on the annihilation cross section to 
a neutrino-antineutrino pair. Comparing this with lim- 
its on the same channel from Ref. [43j], whose technique 
we follow, suggests that the limits presented in this work 
would not be substantially strengthened by the SuperK 
data set, certainly not to the point where the neutrino 
constraints on (u(T)Bicm were competitive with any of the 
stronger bounds. 



VI. DISCUSSION 

Fig. |4] collects our upper limits on the bremsstrahlung 
cross section (t'(T)Brcm, as calculated in the previous 
sections. The constraint from the antiproton ratio is 
stronger than that from the positron data by a factor of 
^ 5. Nature provides a unique value for {va). Therefore, 
if the bremsstrahlung process saturates the allowed an- 
tiproton limit, then the same process produces positrons 
at a rate down from the observed value by about a fac- 
tor of five0 Gonversely, if the observed positron fraction 
were attributed to the bremsstrahlung process, then the 
same process would overproduce antiprotons by about a 
factor of five. 

It is important to note that the observed antiproton 
flux and ratio are well reproduced by standard astro- 
physical processes, leaving little room for a DM contri- 
bution0 We have not attempted to model this standard 
background, so constraints from antiprotons are likely to 
be significantly stronger than presented here. 

Annihilation to or t'^t~ is not as helicity- 

suppressed as to electrons. Even so, the helicity 
factors which suppress the s-wave are (m^/M^) ~ 
10-^ X (Mx/300GeV)-2 and [rrir/M^f ~ 3 x 10"^ x 
(M^/300 GeV)^^, which are comparable to the factor by 
which the p- wave is velocity suppressed, ^ 10^^. Since 
bremsstrahlung overcomes both suppressions, it can also 
be important for annihilation to muons and tans. (And, 
of course, any helicity suppression is especially strin- 
gent for annihilation to and v-^ as rrii, ~ 0.) In the 
case of relatively light DM annihilating to taus, the he- 
licity suppression is not as pronounced. Furthermore, 
W / Z bremsstrahlung may be kinematically forbidden. 
In any case, we note that annihilation to r can never be 
purely leptophilic, as the r has significant hadronic decay 
modes. 

Note that in the model we consider, emission of pho- 
tons or massive gauge bosons can always lift helicity sup- 
pression, however the effect will be much greater in mag- 
nitude when the DM and scalar exchange particles are 
nearly degenerate in mass (such as the co-annihilation re- 
gion of mSUGRA) . In the present work we consider such 
a region of parameter space. Specifically, we present re- 
sults for /X = {Mn/M^f = 1.2, though our conclusions 
remain valid for any value of ji where the bremsstrahlung 
processes dominate the 2 to 2 body processes. As can be 
seen in Ref. [l^] , the rates for these two processes become 
comparable at /j, ~ 10. 

Consider now scenarios where DM annihilation to a 
lepton pair is not helicity-suppressed. As examples, one 



^ Note that we have not compared the spectral shapes of the DM 
signals with those of the observed fluxes, nor tried to fit the data. 

^ Reference l5lll notes that a highly anisotropic diffusion model, as 
might be invoked to accommodate galactic winds, may suppress 
the antiproton flux to a value possibly below the PAMELA flux. 
We do not consider anisotropic diffusion in this work. 
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FIG. 4. Upper limits on (i'cr)Brcm using the 'med' diffusion parameter set. Shown are constraints based on the Fermi EGB 
(solid, red), + e~ flux (dots, orange), e^/(e^ + e~) ratio (dashes, green), p flux (dot-dashes, blue), p/p ratio (dot-dot- 
dashes, magenta), H.E.S.S. gamma rays (dot-dot-dot-dashes, maroon), and neutrinos (dot-dash-dashes, cyan). Also shown for 
comparison is the expected cross section for thermal relic dark matter, 3 x 10"^" cm'^ /s (dot-dot-dash-dashes, black). 
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FIG. 5. As for Fig. (4] using the 'min' (left) and 'max' (right) diffusion parameter sets. 



may have Majorana DM annihilating via an exchange 
of a pseudoscalar or scalar (the latter is still velocity- 
suppressed at the DM vertex) or Dirac DM annihilating 
via the exchange of a vector, or one may have scalar DM 
annihilations. In these cases, there will still be a s^ignal 
from electroweak bremsstrahlung emission [21, al- 

though it will no longer be the dominant channel. Even 
so, the W/Z decay products can still lead to restric- 



tive constraints. Reference [21] considered an example 
(exchange of a scalar) where EW bremsstrahlung makes 
only a subdominant contribution to the total DM anni- 
hilation rate. In this model, the main contribution to 
the annihilation rate comes from the 2-body annihila- 
tion channels, thus the monoenergetic and neutrino 
fluxes dwarf the gauge boson fragmentation products. 
Nonetheless, Ref. 27| found that the antiproton data still 



9 



provide the most stringent cross section constraints for 
certain parameters. Note however, that the models of 
[2^ \2T\ exphcitly break gauge invariance. A detailed, 
model- independent, treatment of weak corrections may 
be found in Ref. [28]. 

It should be noted that the results presented here are 
not due to an exhaustive survey of all possible DM pro- 
files and parameters. Uncertainties arise from the vari- 
ous choices made in order to present illustrative results. 
In most cases, we have made conservative choices for 
these parameters such that alterations to these selections 
should strengthen the results. In calculating the flux of 
protons, antiprotons, electrons and positrons, all the as- 
trophysical parameters are encoded into a numerically-fit 
function 47] with propagation parameters which are con- 
sistent with a 'median' flux 1^56]. However, by assuming 
alternate parameters, e.g. from the 'max' or 'min' fiux 
scenarios, our results may be strengthened or weakened 
by up to an order of magnitude, as shown in Fig. [5l Our 
conclusions hold in all cases considered, but for the ex- 
treme choice of 'min' diffusion parameter set, where the 
e'^ /{e'^ -j- e~) limits become comparable to those for p/p 
and the e+ + limits become comparable to those for 
gamma rays. 

Obviously, the choice of profile can have a large ef- 
fect on the parameter ranges, and we have adopted the 
NFW profile with = 0.39 GeV/cm-^ throughout this 
work. If one considers nonspherical profiles or dark discs 
then the uncertainty in the value of the local DM den- 
sity may be expanded to accommodate a value between 
0.2 GeV/cm-3 and 0.7 GeV/cm-^ji^. Note, however, 
that changes to the DM profile would move all the pre- 
dicted fiuxes, and thus the corresponding cross section 
constraints, in the same direction. For the calculation of 
the extragalactic fluxes, the cosmic source clustering fac- 
tor (/(z) of Eg. dT^ ) can vary by an order of magnitude 
depending on the profile and inclusion or exclusion of 
subhalos 4iJ] (for example, in the case of the Moore pro- 
file [52], our choice for normalization could be increased 
by a factor of 10 [1^). This would lead to tighter con- 
straints coming from the gamma-ray signals. 

We have also neglected the signals produced by inverse 
Compton scattering (gamma rays) and synchrotron radi- 
ation (radiowaves) of the electrons and positrons as they 
propagate in the galaxy (see, e.g., Ref. [13, HI]). Note, 
however, that these effects are properly included in the 



electron energy- loss formalism we adopt. We expect the 
gamma rays produced directly from the annihilations to 
dominate the constraints. Additional gamma rays from 
inverse Compton would only strengthen our results (but 
make them less robust). 

VII. CONCLUSIONS 

If DM is Majorana in nature, then its annihilation 
to fermions may be suppressed due to helicity consid- 
erations. However, when the DM mass is greater than 
M\y/2, both electroweak and photon bremsstrahlung 
may lift this suppression, thereby becoming the domi- 
nant channel for DM annihilation. This permits the in- 
direct detection of models for which the annihilation cross 
section would otherwise be too suppressed to be of inter- 
est. Subsequent decay of the emitted W and Z gauge 
boson will produce fiuxes of electrons, positrons, neutri- 
nos, hadrons, and gamma rays. The aim of the present 
work has been to study the spectra of these particles 
as a tool for indirect detection of DM. By comparing 
these fiuxes to cosmic ray data we have been able to con- 
strain the DM annihilation cross section in such models. 
From these constraints, we find that the observational 
data pertaining to the flux of antiprotons combined with 
those of positrons make it difficult for helicity-suppressed 
2 — > 2 leptophilic DM annihilation to be the source of the 
recently detected cosmic ray anomalies. For these mod- 
els, the bremsstrahlung processes dominate. The primary 
culprit is the hadronization of the gauge bosons, which 
leads to a significant antiproton flux. This result high- 
lights the difficulty of producing lepton-only final states 
even in a model expressly designed for just such a pur- 
pose with 2—^2 annihilation. 
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